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Abstract 

On the triangular lattice, for J2/J1 between 1/8 and 1, the classical Heisen- 
berg model with first and second neighbor interactions presents four-sublattice 
ordered ground-states. Spin- wave calculations of Chubukov and Jolicoeur 
and Korshunov |jT^ suggest that quantum fluctuations select amongst these 
states a colinear two-sublattice order. From theoretical requirements, we de- 
velop the full symmetry analysis of the low lying levels of the spin-1/2 Hamil- 
tonian in the hypotheses of either a four or a two-sublattice order. We show 
on the exact spectra of periodic samples {N = 12, 16 and 28) how quantum 
fluctuations select the colinear order from the four-sublattice order. 
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I. INTRODUCTION 



Symmetry breaking and the selection of a particular macroscopic state amongst many 
degenerate ones result in part from infinitesimal external causes. In the case of planar Neel 
order, the plane of antiferromagnetic ordering, for example, is chosen by the environment, 
whereas the possibility of antiferromagnetic symmetry breaking and the nature of the anti- 
ferromagnetic order are intrinsic and deeply rooted in the spectral properties of the low lying 
levels of the Hamiltonian on a given lattice ■ Two features have to be considered in this 
respect: "the ground-state" and the first excitations of the system. In the past, interest has 
mainly been focussed on the "first excitations": the so-called antiferromagnetic magnons. 
The interest in the ground-state has been limited to the measurement of the order parameter 
modulus. The approach of this problem through exact diagonalizations on small samples 
has led us to focus on the nature of this "ground-state": the eigenstates of the Heisenberg 
Hamiltonian on a finite lattice of sites are eigenvectors of total spin S and, in all presently 
studied cases, the absolute ground-state is 5 = or 5* = 1/2 (depending on the number of 
sites in the sample) . If we consider the even site samples, the S = absolute ground-state is 
spherically symmetric: it does not break the rotational symmetry of the Hamiltonian and as 
such is insufficient to describe a Neel antiferromagnetic state. As underlined by Anderson in 
1952 [|I|, the Neel symmetry breaking state arises from a linear combination of a macroscopic 
number of levels {E} with different S values which in the thermodynamic limit collapse to 
the absolute ground-state faster than the softest magnons. 

This set of levels {E} - called QD JS for Quasi Degenerate Joint States in - has spe- 
cific symmetry and dynamical properties which embody the characteristics of the symmetry 
breaking phase. Let us recall that for a finite solid, the low lying levels are eigenstates of the 
total momentum and they indeed collapse to the ground state in the thermodynamic limit 
faster than the softest phonons. This is a wave packet of these eigenstates that localizes the 
center of mass. Here, in an ordered antiferromagnet, the multiplicity {E} is associated with 
the dynamics of the order parameter. In other words, the knowledge of the symmetry and 
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dynamical properties of this set of eigenstates yields the nature of the ordered phase. For 
the case of a 2d-Neel phase, this set of levels {E} collapses to the ground-states as A^~^, that 
is, faster than other quantities, in particular faster than the softest magnons which converge 
to the ground-state as N~^^'^ [^]. Understanding the symmetry and dynamical properties 
of these low lying levels of the Heisenberg Hamiltonian on the triangular lattice leads to a 
consistent picture of an ordered ground-state with three sublattice Neel order; this reconciles 
spin- wave theories and exact diagonalizations approaches |@,^. 

More subtle symmetry breakings still exist when two or more different kinds of order are 
classically degenerate. In the pure classical case. Villain et al have shown that thermal 
fluctuations could select a specific order. The selected order has softer excitation modes and 
therefore, for a given low energy, a larger density of excitations and a larger entropy: Villain 
et al called this mechanism "order by disorder". This concept has been rather fruitful 
for the studies of classical and quantum antiferromagnets [^0. 

The existence of competing interactions is indeed the main cause of classical ground- 
states degeneracy. As a generic example, one can consider the so-called Ji — J2 model on 
a triangular lattice with two competing antiferromagnetic interactions. This Hamiltonian 
reads: 

H = 2Ji ^ Si.Sj + 2J2 ^ Si.Sfc (1) 

<i,j> «i,k» 

where Ji and J2 = aJi are positive and the first and second sums run on the first and second 
neighbors, respectively. The classical study of this model has been developed by Jolicoeur 
et al 0. They have shown that for small a (a < 1/8) the ground-state corresponds to 
a three-sublattice Neel order with magnetizations at 120° from each other, whereas for 
1/8 < a < 1, there is a degeneracy between a two-sublattice Neel and a four-sublattice 
Neel order (see Fig.|l]). Chubukov and Jolicoeur and Korshunov |T^ have then shown 



that quantum fluctuations (evaluated in a spin-wave approach) could, like thermal ones, lift 
this degeneracy of the classical ground-states and lead to a selection of the colinear state 
(see Fig.|l] ). As usual for spin-1/2 systems, the validity of the spin-wave theory has to be 



checked. The first study of the exact spectrum of Eq.[l| done by Johcoeur et al was not 
incompatible with this conclusion, but was insufficient to yield it immediately. Deutscher 
and Everts found good agreement between spin-wave results for the colinear state and 



exact diagonalizations but their sample geometries were too restricted to fully accommodate 
the four-sublattice order. We show in this paper that a study of the complete dynamical 
"ground-state multiplicity" leads to this conclusion. 

In order to understand the origin of this thermodynamical multiplicity we first study 
exactly solvable models which display either four-sublattice order or colinear order (section 
II). Then, on exact spectra of small samples, we show how quantum fluctuations of increasing 
wavelength select the colinear order (section III). 

II. EXACT SOLVABLE QUANTUM MODELS OF ORDERED SYSTEMS 

These model Hamiltonians are obtained by retaining the Fourier components of the 
Heisenberg Hamiltonian which are compatible either with the two or the four-sublattice 
order. In Fourier components, the Heisenberg Hamiltonian (Eq.|T]) reads: 

Eva 1 
Sk-S_k 7k + — (cosk.(2ui + U2) + cosk.(ui + 2U2) + cosk.(u2 - ui)) , (2) 

k ^ ^ 

where Sk = ^-Sjexpik.rj and 7k = l/3^^cosk.u^ (u^ are three vectors at 120° from 

each other, connecting a given site to first neighbors). In Eq.|^, the k-components associated 

with the k-vectors which keep the sublattices invariant provide the essential features of the 

dynamics of the order parameter. We successively study the case with four-sublattice order 

and the case with two-sublattice colinear order. 

The four vectors which keep the four-sublattice order invariant are k = and the three 

middles of the Brillouin zone boundaries (called in the following k/, liH and k^). In this 

study, we will exclusively consider finite samples with N = 4p sites and with periodic 

conditions: these samples do not frustrate the four (nor the two) sublattice order and they 

effectively present the above-mentioned k vectors in their Brillouin zone. It is straightforward 

to write the contribution of these Fourier components to Ti in the form: 



"Ho = ^(Ji + J2) (S^ - - S| - - Sy , (3) 

where S is the total spin operator and the are the total spin operator of each sublattice. 
^Hq, S^, S^, S^, and S|, form a set of commuting observables. The eigenstates of ^Hq 
have the following energies: 

'E{S, Sa, Sb, Sc, Sd) = f ( Ji + J2) [S{S + 1) - Sa{Sa + 1) 

-SsiSs + 1) - SciSc + 1) - SoiSo + 1)] (4) 

where the quantum numbers Sa, Sb, Sc, Sd run from to and the total spin results 
from a coupling of the four spins 5*^, 5*^, Sc, Sd. 

The low lying levels of Eq.| are obtained for Sa = Sb = Sc = Sd = N/S: 

'Eo{S) = -^^^(iV + 8) + ^(Ji + J,)S{S + 1). (5) 

These states, which have maximal sublattice magnetizations S\ = S]^ = Sq = S^) = 
{N/8 + l)A^/8, are the rotationally invariant projections of the bare Neel states with four 
suhlattices. This is the single physical origin of all properties of {'^-E'}. These levels have 
an energy collapsing to the absolute ground-state as iV~^ justifying the name of tower of 
states or "ground-state multiplicity" given to {^E}. In this exactly solvable model there are 
no quantum fluctuations to renormalize the sublattice magnetization; quantum fluctuations 
will be introduced by the discarded part of the Hamiltonian (Eq.^. 

As we will now show, this multiplicity {^E} can be entirely and uniquely described by 
its symmetry properties under spin rotation and transformation of the space group of the 
lattice. 

Let us begin by the SU(2) properties induced by the fact that these states represent 
the coupling of four identical spins. The degeneracy of each 5* level is (25* + l)Ns where 
the factor (25* + 1) comes from the magnetic degeneracy and Ns is the number of different 
couplings of four spins, each of length A^/8, giving a total spin S. This number is readily 
evaluated by using the decomposition of the product of four spins 72/8 representations of 

su{2) {v^/^y. 
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{^E} = ® ® P^/s ® (6) 

in spin S irreducible representations (T)^). One obtains: 

Ns = l ( -35^ + S{N + 1) + 2 + ^ ) for S <^ 
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(7) 

Kt— ) 

Note that this degeneracy depends both on 5* and and not only on the total spin S 
as is the case for two or three-sublattice problems. In fact, in the latter two cases, which 
describe Neel order on a square or triangular lattice, the objects to be considered stem 
from the coupling of two or three angular momenta: they have perfect counterparts in the 
orbital three dimensional world which are rigid rotators and tops with well known quantum 
numbers, depending only on S. More generally, a Neel order with p sublattices on a finite 
sample of spins gives rise to a "ground-state multiplicity" of the order of A^^. 

The determination of the space symmetries of these eigenstates allows a complete speci- 
fication of {'^E}. The four-sublattice order is invariant in a two-fold rotation {71-,^): thus the 
eigenstates of {'^E} belong to the trivial representation of C2. As it arises from the coupling 
of four identical spins, this subset of levels forms a representation space of 5*4, the permuta- 
tion group of four elements. The eigenstates of {'^E} could thus be labeled by the irreducible 
representations of 5*4 (see Table |). Indeed, the complete analysis of all the eigenstates of 
Eq.|l| is usually done through the more general point of view of the space group of the lat- 
tice. But it is straightforward to show that in the four-sublattice subset of solutions, each 
element of the space group maps onto a permutation of 5*4: one step translations map onto 
products of transpositions as {A, B){C, D), three-fold rotations onto circular permutations 
of three sublattices {A, B, C) and so on. The complete mapping of the space symmetries 
of the four-sublattice order onto the permutations of 6*4 is given in Table | together with 
the character table of S^. Each irreducible representation of 6*4 can thus be characterized 
in terms of its space symmetry properties. As noted above they are all invariant in 71.-,^. 
Analysis of the properties under translation shows that ri,r2 and correspond to the 



wave-vector 0, whereas r4 and Fs have a wave-vector k//, k/ or kj. Fi and F2 are invariant 
under C3, whereas F3 is associated with the two complex representations of this same group. 
Finally, Fi and F4 are even under axial symmetry whereas F2 and F5 are odd. The number 
of replicas of F, that should appear for each 5* is then computed in the S, Ms subspace with 
the help of the trace of the permutations of S4: 



n 



(S) 



^J2TT{R,\s)Mk)N,i{k) 



(8) 



where Rk is an element of the class k of S4, Nc\{k) is the number of elements in this class 
and Xi{k) is the character of the class k in the irreducible representation Fj (see Table |). 
The values of the traces for a given total spin S are then found as: 



Tl^((AB.C)P=T.((A,S.C)[^^J-T.((AB 



C) 



Ms=S+l 

In each Ms subspace of {'^E}, it is straightforward to find the trace of the elements of S4: 
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Tr /, 



Ms, 

Tt[{A,B){C,D) 



Ms 



Ms 



TTi^{A,B,C) 
Tr ({A,B 



Ms 

Tr ( {A, B, C, D 



Ms 



N/8 

^t+v+x+y,Ms 

t,v,x,y=-N/8 
N/8 

v,Ms 

t,v=-N/8 
N/8 

^3t+v,Ms 

t,v=-N/8 
N/8 

^2t+v+x,Ms 

t,v,x=-N/8 
N/8 

t=~N/8 
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where t,v,x,y are the 2;-components of the total spin of each sublattice (constrained to 
vary between and —N/8) and 6ij denotes the Kronecker symbol. Using Eqs.p|j9|JT0| 
one readily obtains the number of occurrences of each Fj for each 5* subset of {^-E} (Table 
|T|). We have thus obtained the complete determination (all quantum numbers, and all the 
degeneracies) of the family of low lying levels describing the ground-state multiplicity of the 
four-sublattice Neel solutions. 



Let us now consider the colinear solutions (Fig.p. They are particular solutions of the 
four-sublattice case and we will rapidly go through the same scheme of analysis, indicating 
mainly the new points. The two vectors which keep the two sublattices invariant are 
and the middle of one side of the Brillouin zone (the vectors k/, hn and k^ correspond 
respectively to the colinear solutions (a), (6) and (c) in Fig.^. Extracting a specific set of 
two wave-vectors from Eq.|^, we find the following contribution to the total Hamiltonian: 



N 

The corresponding low energy spectrum for Sa = Sf3 = N/4 is: 



(11) 



'Eo{S) = -^^^{N + 8) + ^{J, + J,)SiS + 1) (12) 

and is degenerate with the four-sublattice low energy spectrum (see Eq.^. Here, the two 
sublattices have maximal spins Sa = Sjj = N/4. These new solutions arise from the sym- 
metric coupling of the spins of two sublattices of the four-sublattice order: Sq = Sa + Sb or 
Sa = Sa + Sc or S^ = Sa + Sd with the counterparts for Sp. As there are three ways to do 
this coupling, the colinear solutions have a degeneracy. The representation space is thus 
the sum of three products V^^^^V^^'^ . It is not a direct sum since V^/'^^A, B) 0V^^^(C, D) 
and D^/^(A, C) ®T>^/^{B^D) have in common the same (symmetric) irreducible represen- 
tation with a total spin N/2. On a A^-sample, the representation space of the ground-state 
of the colinear solution is: 

The degeneracy is thus 3(25* + 1) for all S values except for S = N/2 where it is {2S + 1). 

As for the four-sublattice order, the space-group analysis is done as for the two-sublattice 
order, but the number of occurrences of each irreducible representations Fj is now different 
since the space {^E} is smaller than {'^E}. The calculation could be done along the same 
lines as for the four-sublattice order. The problem, however, is much simpler because for 
each S value there are only three replicas of arising from the Z3 group (Eq.^ and Fig.p. 
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This allows direct computation of the permutation traces in each S subset of {'^E}. Using 
the coupling rules of two angular momenta (and in particular the fact that the S eigenstate 
resulting from the coupling of two integer spins changes sign as (—1)'^ with the interchange 
of the two parent spins) one obtains (for S ^ N/2): 



Tr 
Tr 
Tr 
Tr 
Tr 



{A,B){C,D) 
{A,B,C) 
{A,B) 



3 
1 

1 



2(-l) 



(14) 



iA,B,C,D) 



Therefore, the colinear solution is simply characterized by Fi and for even S and r4 for 
odd S. 

From these equations (Eqs.p|,P JTO| , p!^ ) , the symmetries of all states of the tower are fully 
determined both for the four-sublattice order {'^E} and for the colinear order {^E}. 

Going back now to the original (Ji — J2) model, we have to account for quantum fluc- 
tuations generated by the discarded part of Ti. This perturbation does not commute with 
sublattice total spins and consequently reduces the sublattice magnetization. Nevertheless, 
it preserves all the symmetries of the Neel state and thus also the ones of the levels of {^E} 
or {^E}. Then the question is : do quantum fluctuations conserve qualitatively the dynam- 
ics of these levels or not ? If these levels remain the low lying ones of the exact spectra 
with overall dynamics qualitatively similar to that of the bare Neel state (Eqs. ^p!2D, then 
the quantum model will be ordered at T = 0. By qualitatively, we mean that the leading 
term of the energy of the exact subset {E} behaves as Pj^{Ji + J2)S{S + 1), where f3 is 
a renormalization factor. This factor is related to the spherical homogeneous susceptibility 
of the sample ||^, even if, in general, the tensor of susceptibilities is not spherical because 
quantum fluctuations lift the degeneracies of {'^E} of the exactly solvable model. 
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III. EXACT SPECTRA OF SMALL PERIODIC SAMPLES 



We have determined the low (and high) energy levels of the Ji — J2 Hamiltonian in 
each irreducible representation of SU{2) and of the space group of the triangular lattice 
for small periodic samples with N = 12, 16 and 28. The spectra are displayed in Fig.|^ 
and Fig.^. We directly see in the upper parts of these figures the set of QDJS ("ground- 
state multiplicity" ) well separated from the set of levels corresponding to the one magnon 
excitations. We have verified that the QDJS form a set of levels with the exact properties 
of the above defined {^E} subset. The action of quantum fluctuations could then be read 
in the lower parts of the figures. As expected, quantum fluctuations lift the degeneracies 
which are present in the exactly solvable model and stabilize the eigenstates with the lower 
S values. Nevertheless, the low lying energies per site still group around a line of equation 
+ 8/35(5 + l)(Ji + J2)/N^ with P = 1.004 (resp. 1.055) for = 16 (resp. 28). The 
number and space symmetries of these levels for each S and N value are exactly those 
required by the above analysis of the four-sublattice Neel order. Moreover, it is already 
visible on the = 16 sample and quite clear on the = 28 sample that a dichotomy 
appears in this family (see Fig.^). The lowest levels of this tower of states appear to be 
Fi, F3 or F4 representations depending on the parity of the total spin. They precisely build 
the family {"^E} of isotropic projections of the colinear solutions given above (Eq.^). 

We see in Fig.^ that the difference between the energy per bond of the colinear states 
and that of the other states of the four-sublattice order roughly increases by a factor 4 from 
the = 16 to the A^ = 28 sample. This strongly suggests that the four-sublattice order will 
disappear in the thermodynamic limit and only the colinear order will persist. This result 
supports the conclusion of the spin-wave expansion ||8|-[T0| concerning the selection of the 
colinear state in the Ji — J2 model for 1/8 < a < 1 for spins 1/2. 
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IV. CONCLUDING REMARKS 



It appears from the two situations that we have studied (triangular Heisenberg model 
and this model) that the symmetry and dynamical analysis of the low lying levels of 
a Hamiltonian likely to exhibit ordered solutions give rather straightforward answer to the 
kind of order to be expected. The method is rapid, powerful and unbiased: it does not 
require any a priori symmetry breaking choice: if a specific order is selected, one should see 
it directly on the exact spectra. Moreover, as it is essentially exact, there are no questions 
relative to the convergence of the expansion as in the spin-wave approach. On the other 
hand, as the sizes amenable to computation are limited, there is, in the exact approach, a 
cut-off of the long wavelength fluctuations. Results so obtained should thus be examined 
in light of a finite size scaling analysis. The present work nevertheless shows that it is not 
necessary to invoke quantum fluctuations with very long wavelengths to select the colinear 
order. 
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TABLES 



5*4 


I 


{A,B){C,D) 


{A,B,C) 


{A,B) 


{A,B,C,D) 


G 


I 


t 


^27r/3 


a 




Nel 


1 


3 


8 


6 


6 




1 


1 


1 


1 


1 




1 


1 


1 


-1 


-1 




2 


2 


-1 








r4 


3 


-1 





1 


-1 




3 


-1 





-1 


1 



TABLE I. Character table of the permutation group S'4. First hne indicates classes of per- 
mutations. Second line gives an element of the space symmetry class corresponding to the class 
of permutation. These space symmetries are: t the one step translation (A ^ B), 'Ri-Kjz (resp. 
7^27r/3) t'^^ three-fold rotation around a site of the D (resp. ii?)-sublattice, and a the axial symmetry 
keeping invariant C and D. N^i is the number of elements in each class. 
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AT = 16 

5 





1 


2 


3 


4 


5 


6 


7 


8 


Tlr^ (S) 


1 





2 





2 


1 


1 





1 










1 






















2 





2 


1 


2 





1 













2 


2 


3 


2 


2 


1 


1 










2 


1 


2 


1 


1 












= 28 
S 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


14 


nrAS)+nr^{S) 


2 





5 


1 


5 


3 


4 


2 


4 


1 


2 


1 


1 





1 




3 





4 


2 


5 


2 


5 


2 


3 


1 


2 





1 








nr^iS) + nr^{S) 





7 


6 


11 


9 


12 


9 


10 


6 


6 


3 


3 


1 


1 






TABLE II. Number of occurrences nr^ (S) of each irreducible representation with respect to 
the total spin S. For N = 28, and as well as nr^ and npg have been added because this 
sample does not present any axial symmetry. 
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FIGURES 

FIG. 1. Top: four-sublattice classical ground-state. Spins in the sublattices A and B, as well 
as spins in C and D, make an angle 26. The plane of the spins of A and B makes an angle (p 
with the plane of the spins of C and D. Bottom: the colinear solutions with the three possible 
arrangements (in this case, classical spins in sublattices A and B are antiparallel) . 

FIG. 2. Top: complete spectrum of the = 16 periodic sample with respect to for 
J2/J1 = 0.7. Bottom: enlargement of the difference between the exact spectrum and the energy 
of the low lying levels of the model Hamiltonian (Eq.|5| or Eq.l2). The ground-state multiplicity 
is well separated from the magnons. 

FIG. 3. Partial spectrum of the = 28 periodic sample (same legend as for Fig.||). Bottom: 
the tower of states of the four-sublattice order {^-E} lays under the dashed line. Above appear the 
first magnons. Above the dotted line are represented the first excited homogeneous states. In the 
magnon multiplicity (k 7^ 0, k//,k/ or kj), for 5 < 5, only the 5 lowest states of each irreducible 
representation have been computed. 

FIG. 4. Enlargement of the A^ = 16 and A^ = 28 QDJS. A global contribution l3Eo{S) is 
subtracted from the exact spectrum. This contribution describes the overall dynamics of the order 
parameter in this finite sample, (3 measures the renormalization of this dynamics by quantum 
fluctuations (see text and Ref. [^). The bars represent eigenstates which belong both to {^E} and 
The triangles indicate states which belong to {^-E} but not to {^E}. We see that, with 
increasing sizes, the tower of states of the colinear order separates from the four-sublattice order. 
For A^ = 28, the two states of {^E} with even S are quasi degenerate and cannot be distinguished 
at the scale of the figure. 
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